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OPTIMAL TIME DECAY OF NAVIER-STOKES EQUATIONS 
WITH LOW REGULARITY INITIAL DATA 

JUNXIONG JIA AND JIGEN PENG 


Abstract. In this paper, we study the optimal time decay rate of isentropic 
Navier-Stokes equations under the low regularity assumptions about initial 
data. In the previous works about optimal time decay rate, the initial data 
need to be small in H^- N ^ 2 ^ 2 (R N ). Our work combined negative Besov space 
estimates and the conventional energy estimates in Besov space framework 
which is developed by R. Danchin. Though our methods, we can get optimal 
time decay rate with initial data just small in B N / 2 - 1 ’ N / 2 + 1 r\B N / 2 - 1 ’ N / 2 and 
belong to some negative Besov space(need not to be small). Finally, combining 
the recent results in m with our methods, we can only need the initial data 
to be small in homogeneous Besov space B N / 2 — 2 > N / 2 n B N / 2 ~ 1 to get the 
optimal time decay rate in space L 2 . 


1. Introduction and main results 

In this paper, we consider the large time behavior of global strong solutions to 
the initial value problem for the compressible Navier-Stokes system 

{ d t p + div(pu) = 0, 

dt(pu) + di v(pu (g> u) + VP(p) — p\7u — (A + /i)Vdi vu = 0, 

(p,u)\ t =o = ( Po,u 0 ). 

The variables are the density p > 0, the velocity u. Furthermore, P = P{p) is the 
pressure function. The viscosity coefficient satisfy p > 0 and A + 2p > 0. p is a 
constant and, for simplicity, we can take p = 1 in the present paper. 

There are many results on the problem of long time behavior of global solutions 
to the compressible Navier-Stokes equations. For multi-dimensional Navier-Stokes 
equations, the H s global existence and time-decay rate of strong solutions are 
obtained in whole space first by Matsumura and Nishida mm and the optimal 
L p (p > 2) decay rate is established by Ponce 0. To conclude, the optimal L 2 time 
decay rate in three dimension is established as 

(1.2) j|(p — p,u)(t)|| i 2( R 3) < C(1 + t) 4 

with (p, 0) to be the constant state and the initial data to be a small perturbation in 
Sobolev space P[ l with l > 3. In [4], the author developed a general energy method 
for proving the optimal time decay rates of the solutions. They removed the small 
requirement about the initial data in L^R 3 ) norm, however, they need the initial 
data belong to the negative sobolev space H~ s with s G [0,3/2). Their method is 
so interesting that it needs not to analyze the Green’s matrix of linearized system. 
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H. Li and T. Zhang 0 introduced the negative Besov space for the initial data. 
More specifically, they need 

II (Po — P) u o)||^fi(R3)ns-^ o (K 3 ) ^ — 4, s £ [0,1] 

small enough. Then they can obtain the faster time decay rate than the classical 
results m as follows 

(T3) ||(p-p,u)(t)|| L 2 (R 3) < (7(1 

Time decay in the energy space imply the dissipative properties of the compress¬ 
ible viscous fluids, however do not describe the information of wave propagation. 
To understand the wave propagation of compressible viscous in multi-dimension, 
Hoff-Zumbrun mw consider the Green’s function of the compressible isentropic 
Navier-Stokes equations with the artificial viscosity and show the convergence of 
global solution to diffusive waves. Liu-Wang [8] investigate carefully the properties 
of Green’s function for isentropic Navier-Stokes system and show an interesting 
convergence of global solution to the diffusive waves with the optimal time-decay 
rate in odd dimension where the important phenomena of the weaker Huygens’ 
principle is also justified due to the dispersion effects of compressible viscous fluids 
in multidimensional odd space. This is generalized to the full system later in [5] 
where the wave motions of other types are also introduced. Therein, the optimal 
L°° time-decay rate in three dimensions is 

(1.4) ||(p(t) — p, u(t))||ioo( R 3) < (7(1 + t) 3//2 . 

The same decay property also proved in the exterior domain maim, the half space 

[12 El- 

Inspired by the previous works, we also introduced the negative Besov space. 
Compared to the above work, our work highly reduced the requirements about 
the regularity of the initial data and through our proof we revealed an important 
difference about critical Besov space framework between compressible field and 
incompressible field. Below, we give our results specifically. 

Theorem 1.1. Assume that (po — p, uo) £ qN/ 2 -i,n/ 2 +i x ^n/ 2 -i,n/ 2 ^n j Qr an 
integer N > 2. Then there exists a constant ao > 0 such that if 

(1.5) ||P0 — p||gJV/2-l,JV/2 + l + ||'Uo|lBW/2-l,iV/2 < Op) 

then the problem dp admits a unique global solution (p, u) satisfying that for all 
t > 0, 

^ Up — pIIc(r+.b n / 2 -un/2+i) + ||m|| C ( R+i bn/2-i,n/2) 

+ Up ~ p\\L' i -(K.+,B N / 2 + i ) + IMIl i (r+,b n / 2 + i ’- w /2+2) — Ca.Q. 

If further, (p 0 — p, u 0 ) £ for some s £ [0, JV/2], then for all t> 0, 

(1.7) \\p{t) + IMOIIs-^ < c 

and 

minU,N/2-l) + s 

(1.8) \\p(t) - p\\g e + \\u(t)\\g e <C(l+t) 2 for — s < £ < N/2. 
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Remark 1.2. For N = 3 and s = |, taking l = 0 in (11.81) . we will have 

II P(t) -p|U= + IK*)IU= < ||p(t) ~p\\bo + l|w(i)llso 

< <7(1+4)"*, 

which coincide with the classical results m- However, we need not the L 1 , H 3 


norm of initial data to be small. 

Remark 1.3. Taking iV = 3, s = § and £ = |, we can obtain the L°°-time decay 
for 3-D problem as follows 


II p(t) ~ P|U°° + ||w(f)||L~ < C\\p(t) — p||^l/2,5/2 + \\u(t)\\g 1/2 ,3/2 

< Cil + t)- 1 . 


This time decay rate is a little slower than the optimal L°°-time decay rate. How 
to get optimal L°°-time decay rate under the low regularity assumptions on initial 
data is our future aim. 

Remark 1.4. The methods we used here seems can be generalized to a lot of other 
fluid dynamic models such as Navier-Stokes Poisson, magnetohydrodynamic flows. 
This will be left to the future work. 

Remark 1.5. Recently, C. Wang, W. Wang and Z. Zhang pi] construct global so¬ 
lution for Navier-Stokes equations allow the initial density and velocity both have 
large oscillation. Combining their results with our proof methods, we can get the 
following results. 

Theorem 1.6. Suppose the equilibrium state p > 0 and dimension N = 3. Let Co to 
be a small constant, P+ = sup Co/4 < p < 4c -i |P (fc) (p)l and P- = inf Co/4 < p < 4c -i \P'{p)\- 
Assume that the initial data (po,uo) satisfies 



where s > 3. There exist two constant ci = ci(A, p, Co, p, P+, P-) and C 2 = 
C 2 (A, p, co, p, P+, P -) such that if 



\\P0 ~ PIIrI/ 2,3/2 < Cl, |u 0 |Ll/2 < 

D 2,l ^ 2,1 


ll u o||/j-s < 


(1 + IlPo - Plls3/2)(1 + ||po - p ||^2 + ||«o||^ 3 ) 


for some 6 G (—A, |), then there exists a unique global solution ( p,u ) to the 1 11 . If) 
satisfying 


P>^j, P ~ P £ C([0, +oo); H s ), 
u G C([0, +oo); II s - 1 ) n L 2 (0, T; H s ), for any T > 0. 


In addition, if c\, C 2 small enough, (po — p, uq) £ 8 2 * f or some s G [0, |], then 
for all t > 0, we have 


\\p(t) - pWb-^ + IKOIIb-^ < C- 


and 


IIpW ~ p\\b* 1 + II u (0II.b* 1 ^ c(i + *) 2 i 
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for — s < £ < where C > 0 is a positive constant depending on the initial data. 

The proof of this Theorem is similar to the proof of Theorem o so we only 
give Remark 13.21 Remark 14.21 and Remark 15.21 in the following to indicate how to 
get this result. 

The present paper is structured as follows. In section 2, we will give some basic 
knowledge about homogeneous Besov space. Then in section 3, we get a differen¬ 
tial type inequality about a hyperbolic parabolic system comes from compressible 
Navier-stokes equations. In section 4, we proved that Navier-Stokes equations con¬ 
serve the navigate Besov norms. Theorem ll.il will be proved by using interpolations 
of Beosv space in section 5. 

Notation. Throughout the paper, c and C stands for a “harmless” constant, 
and we sometimes use the notation A < B as an equivalent of A < CB. The 
notation A « B means that A < B and B < A. 


2. Preliminaries 

In this section, we will give some basic knowledge about Besov space, which can 
be found in 1141 . The homogeneous Littlewood-Paley decomposition relies upon a 
dyadic partition of unity. We can use for instance any </> £ supported in 

C := {£ £ R^, 3/4 < |£| < 8/3} such that 

5>(2" fl 0 = l if 

? ez 

Denote h = we then define the dyadic blocks as follows 

A q u := (f{2~ q D)u = 2 qN f h(2 q y)u(x — y) dy, and S q u = Afcit. 

The formal decomposition 

u = A q u 


is called homogeneous Littlewood-Paley decomposition. The above dyadic decom¬ 
position has nice properties of quasi-orthogonality: with our choice of <f>, we have 

AfcA q u = 0 if |A; — q\ > 2, 

A k (Sq-iA q u) = 0 if \k — q\> 5. 

Let us now introduce the homogeneous Besov space. 


Definition 2.1. We denote by S' h the space of tempered distributions u such that 

lim SjU = 0 in S'. 

j-*-oo 


Definition 2.2. Let s be a real number and (p, r) be in [1, oo] 2 . The homogeneous 
Besov space consists of distributions u in S' h such that 



V 


r 


J2^ s \\Afu\\ r LP 


< Too. 
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From now on, we denote B£ 1 as B * and the notation B s will stand for i3| 1 . 
The study of non stationary PDE’s usually requires spaces of type U T (X) ■= 
L r (0,T-,X) for appropriate Banach spaces X. In our case, we expect X to be a 
Besov spaces, so that it is natural to localize the equations through Littlewood- 
Paley decomposition. We then get estimates for each dyadic block and perform 
integration in time. However, in doing so, we obtain bounds in spaces which are 
not of type L r { 0, T ; B s ). This approach was initiated in .19] and naturally leads to 
the following definitions: 

Definition 2.3. Let ( r,p) £ [l,+oo] 2 , T £ (0,+oo] and s £ R. We set 

:= J2 2qs 

and 

L r T (B°) := {u £ V t {B% < +^} . 

Owing to Minkowski inequality, we have L r T (B s ) LIp(B s ). That embedding is 
strict in general if r > 1. We will denote by Ct(B s ) the set of function u belonging 
to L™(B S ) nc([o,r] ; B s ). 

Another important space is the hybrid Besov space, we give the definitions and 
collect some properties. 



i /r 


Definition 2.4. Let s,t £ R. We set 


and 


||u||:= £ 2 qs \\A q u\\ LP + Y, °- qt II A,«||lp. 

q< 0 q>0 

B°/(R n ) := {u £ S' h (R N ), \\u\\ h s,t < +oo| . 


Lemma 2.5. 1) We have Bp’ s = B 

2) If s <t then B*’* = 5*05*. Otherwise, Bp’* = Bf, + Bp. 

3) The space Bp ,s coincide with the usual inhomogeneous Besov space. 

4) If s l < S 2 and ti > t 2 then Bp 1 ’* 1 Bp 2 ’* 2 . 

Throughout the paper, we shall use some paradifferential calculus. It is a nice 
way to define a generalized product between distributions which is continuous in 
functional spaces where the usual product does not make sense. The paraproduct 
between u and v is defined by 

T u v := Y, S q -iuA q v. 

We thus have the following formal decomposition: 

uv = T u v + T v u + R{u, v), 


with 

R(u , v) := Y^ A q uA q V, A q \= A g _l + Aq + A ?+i . 


We will sometimes use the notation Tfv := T u v + R(u, v). 
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In the last part of this section, we collect some useful lemmas which will be used 
in the sequel. 


Lemma 2.6. [II] A constant C exists which satisfies the following properties. If 
S\ and S 2 are real numbers such that si < S 2 and 6 £ (0,1), then we have, for any 
(p , r) £ [1, oo] 2 and any u £ S' h , 


1^11 D es i+( i_e ) s 2 — II^IIr s i IMU*2 > 

±->7} r Up -p Ur>.r 


Ml • e ai+ (i_0) S2 < 


c 


S 2 — s\ \9 1 — 0 


il -9 


Lemma 2.7. [ 11 j There exists a constant C such that for any real number s and 
any ( p,r ) £ [l,oo] 2 , we have, for any (u,v) £ L°° x B pr , 

II^«' u II_b» r — CIMUMMIb*,.- 

Moreover, for any ( s,t) £ lx (—oo,0) and any (p,ri,r 2 ) £ [l,oo] 3 we have, for any 
( u,v ) £ B ^ )ri x B s pr2 , 

\\ t uv\\b°p < j\\u\\ 

Remark 2.8. Using similar methods as in the proof of Lemma |2.71 we can get the 
following estimates 


Ml n s with - = min < 1 

D p,r 2 f 


1 

ri 


1 

T2 


Lemma 2.9. p3] A constant C exists which satisfies the following inequalities. Let 
(si,S 2 ) be in R 2 and (pi,P 2 Ui; r 2 ) be in [l,oo] 4 . Assume that 

1 1 1 ,1 1 1 

- =-1-<1 and - = -I-<1. 

p pi P2 r ri r2 

If si + S 2 is positive, then we have, for any (u,v) in B p \ ri x B p 2 r2 , 




C 


si+s 


Sl + S 2 


IM 


R s 2 
±J P2 > r 2 


When r = 1 and si + S 2 > 0, we have, for any ( u,v ) in B p 1 i ri x B p 2 r2 , 


||i?(u,u)|| 


s l + s 2 


<Cl|t 


II B “ 


Lemma 2.10. [TB] Let F be an homogeneous smooth function of degree to. Suppose 
that < si,ti, S 2 ,t 2 < 1 + ^. The following estimate hold: 


\(F(D)A q {u-X7a)\F(D)A q a)\ 

< a q 2 ^ ^||u||^jv/2+i \\a\\ Bsi , a2 \\F(D)A q a\\ L2 , 

with '}2 qe %ot q < 1 and 

(f a ’P{q) = a if q < 0, 

r’ fl (q)=P if g > 1 - 
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3. Estimates about hyperbolic parabolic system 


This part is inspired by R. Danchin’s seminal work |15j . however, we need a 
differential type estimates different to the integral type estimates proved by R. 
Danchin. Firstly, let us give some equivalent forms of the systems ED- Denoting 
a = p — 1, supposing P'{ 1) = 1, N > 2 stands for dimension, then we have 


(3.1) 


dta + u • Va + diva = —adivu, 

d t u + u • Vu — pAu — (A + p)Vdivu + Va 

= — h(a ) (pAu + (A + p)Vdivit) — /(a)Va, 
(a,n)|t=o = (po - l,uo), 


where 


M«) = /(a) 


P'jl + a) 

1 -|- CL 


P\ !)• 


Denoting d = A 1 divn and D = A 1 curln with (curfo)l = djZ 1 — diZ-’ , we will 
obtain the following equivalent form 


{ d t a + u • Va + A d = F, 
d t d + u • Vd — vAd — Aa = G, 
d t fl — p£l = H , 
u = — A~ : Vd — A -1 divfl, 


where v = A + 2p and 


F = —adivu, 

G = u ■ Vd + A _1 div (—a • Vu — h(a)(pAu + (A + pjVdivzi) — /(a)Va), 

H = A~ 1 curl (—u ■ Va — h(a)(pAu + (A + p) Vdivu)). 

The equivalence is not so obvious, considering the proof can be found in mum, 
we will omit it. Denote 

E s =L 1 (R+;B S x R s+1 > s+2 ) n C , fe (K + ; B s “ 1>s+1 x B s ~ 1,s ). 

Then, we give the main results in this section. 


Proposition 3.1. Let (a, d, Q) be solution of (13.211 on [0, T), 2 — N/2 < s < l+N/2, 
V(t) = fg ||'a( , r)||^ i v/ 2 + i dr, and V(oo) is bounded. Then there exists a functional 

(3.3) £ s (a,d,f2) ~ ||a(£)||_g a — 2 ,, + ||d(t)||^ a _ 2 , s -i + ||D(t)||^ 3 _ 2 , s -i, 


such that 

d_ 

, dt 
(3.4) 


( e~ cv ^£ s (a , d, fi)) + ce~ cv W (j|a(f)||^ + ||d(i)[|^„ +1 
+ ||fi(t)||*... +1 ) < Ce~ cv W (||F(t )|| 5s _ 2 , 3 + ||(G(t), id(t))||^ s _ 2 , s _ 1 ), 


where C depends on p, A, N , s. 


Proof. Let I\ > 0 to be a constant which will be specified later, define 

(a, d, D) = e~ KV ^ (a, d, fl), 

(F,G) = e~ KVW {F 1 G). 
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Applying the operator A q to the first two equations in (13.21) . we infer that (A g a, A q d) 
satisfies 

df A q a + A q (u • Vd) + AA q d = A q F — KV'(t)A q a, 

dtA q d + A q {u • Vd) — vAA q d — AA q a = AG — KV' (t)A q d. 

Step 1 : low frequencies. Suppose that q < qo and qo is a constant specified in 
Step 2, define 

fq = ll A 9 a||| 2 + ||A g dj|| 2 - r(AA q d\A q d). 

Let oo = | and rao < 1. According to the Holder’s inequality, we have 
|(AA g d|A g d)| < o 0 (||A g a||| 2 + ||Ad||| 2 ) . 

Hence, we know that 


f q * ||A,5||£ 2 + ||A,d||£ a 

Using standard energy methods, we can obtain the following basic energy estimates 

^ 6) 7^ll A 3«ll! 2 + ' Vd)|A g d) + (AA q d\A q a) 

= {A q F\A q a)-KV'(t)\\A q d\\ 2 L ,, 

(3 7) “II VUi, + ■ Vd)|A g d) - v{AA q d\A q d) 

- (AA q a\A q d) = (A q G\A q d) - KV'(i)\\A q d\\ 2 L2 , 


^(AA q a\Ad) + ||AA g d ||| 2 - ||AA g d ||| 2 + (AA,(« ■ Va)|A g d) 
( 3 - 8 ) + (A,(u • Vd)|AA g d) + ^(A 2 A g d|AA g a) = (AA g F|Ad~) 

+ (AG|AA g a) - 2 KV'(t)(AA q a\A q d). 

Performing the following calculation 

CUD + HMD - ^(ESD, 

we have 

l -j t f q + M 1 =M 2 + M 3 -KV'(t)fq, 

where 

( V \ ~ V 7/r* ~ 

" - 2) H AZ Mli 2 + 2 H AA «“I ^ 2 “ — (A 2 A g d|AA g d), 

M2 = (AqFjAqd) + (AqGjAqd) - ^(AA q F\Ad) - ^(A g G|AA g d), 


A4 3 


Since 


-(A q (u ■ Va)|A g a) - (A q (u • Vd)|A g d) + -(AA,(u • Va)|A g d) 
+ ^(A q (u ■ Vd)\AA q d). 


— |(A 2 A g d|AA g d)| < 


,2 „2 


rv a\ 


■||AA,d|| 


2 

L 2 


77 || AA g a| 


2 

L 2 ' 


2 
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Taking r so small such that 

v _r__r^ 

2 2 

Then the term Mi can be bounded from below by c2 2q f q . According to Holder’s 
inequality, we obtain 

M 2 < \\A q F\\ L 2 \\A q d \\ L2 + ||A,G|| La ||A 9 d || L3 

+ r|| AA,F || L 2 1 | A q d\\ L 2 + r|| A 9 G || l2 || AA g a|| L3 
<a q 2 -^- 2 )/ q ||(F,G)|| 5s _ 2 , 3x5fl _ 2 , s _ 1 . 

Using Lemma T2. 1 01 directly, we get the estimates of .M 3 as follows: 

Ms <a q 2- q ^- 2 W\t)[\\a\\^, s + (||A,5|| ia 

+ ||A,d|| La ) +a g U'(t)r( 2 -^- 2 )||AA 9 a|U 2 ||Jb s _ 2 , s _ 1 

+ 2 -^- 2 ) 2 ^||a|| 5 s _ 2 , 3 ||A ? J|U 2 ) 
<a g 2 -^- 2 )U'(i)/ 9 ||(a,J)|| 5s _ 2 , Sx ^- 2 , s - 1 . 

So for the low frequency part, we finally have 

+ c2 2 “f q < Ca q 2-^-V(\\(F, G)|| 5> _ 3 ,. x5 ._ a ,._ 1 
(3.9) dt v 

+ y'(t)||(5,d‘)||B.- 2 ,. X B.- 2 ..- 1 ) -CKV'(t)f q . 

Step 2 : high frequencies. Suppose q > qo + 1 in this part. Define 
/ 2 = ||A 2 A g a ||| 2 + ^||AA g J ||| 2 - ^(A 2 A g a|AA g d). 

Since 

“l(A 2 A g a|AA g d)| < i||A 2 A g h|| 2 2 + ^||AA g J||| 2 , 

we know that 

fq ~ l|A 2 A g a ||| 2 + ||AA g d|| 2 2 . 

Similar to the low frequency part, we have the following basic energy estimates. 

(3 1Q) i|||A 2 A g h||| 2 + (A 2 A g ( M • Va)|A 2 A g a) 

+ (A 3 A g d|A 2 A g a) = (A 2 A g A|A 2 A g a) - A'U'(t)||A 2 A g h|| 2 2 , 

(3 11) 2 dt (AA q (u ■ Vd)jAA q d) + i/||A 2 A g a ||| 2 

- (A 2 A g a|AA g J) = (AA g G|AA g J) - ifU'(t)||AA g d|| 2 2 , 

4(A 2 A g h|AA g J) + ||A 2 A g J || 2 2 + (A 2 A q (u ■ Va)|AA q d) 

+ (AA g (u ■ Vd)|A 2 A g a) + p(A 3 A g J|A 2 A g a) — ||A 2 A g a||| 2 
= (A 2 A g A|AA g d) + (AA g G|A 2 A g a) - 2A'U'(t)(A 2 A g h|AA g J). 


(3.12) 
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Performing the following calculation 


we will obtain 


where 


(MD + 4™ - -<£H, 


l -j t f 2 q + M 1 = M 2 + Ms - KV' 


Mi = i||A 2 A,d||| 2 + i||A 2 A ? 5||| 2 - -J^ 2 A q &\AA q d) 


M 2 =(A 2 A q FjA 2 A q d) + —^(AA q G\AA q d) 

- ^(A 2 A q F\AA q d) - ^(AA q G\A 2 A q d), 


M 3 = - (A 2 A q (u ■ Vd)\A 2 A q d) - A(aa q (u ■ Vd)\AAd) 

+ — ^(A 2 A q (u ■ Va)|AA g d) + (AA q (u ■ Vd)|A 2 A g d)^ . 
We can easily obtain 

|(AA ? d|A 2 A g a)| < ^||A 2 A g h|| 2 2 + i||AA 9 d||l 2 , 

||A 2 A 9 d|| 2 2 >^2 2 (®+ 1 )||AA 9 d||l 2 . 


Taking qo large enough such that 

1 Q ? 

-- 2 2q ° - 4 > 0 . 
v 4 v 3 

Aii can be bounded from below by c/ 2 . Similar to the low frequency part, we have 


M 2 < Ca q f q 2~ q ^ s ~ 2 ^ ||(F, G)|| §s-2,s x §b-2,b-i , 

M 3 < Ca q f q 2- q ^V\t)\\(d, d)||s.- a ..- 1X B.- 2 ..- 1 


So finally, for high frequency part, we know that 
d 


(3.13) 


dt 


f q + cf q <Ca q 2-^~ 2 M\\{F,G)\\B‘- 


2 - s xB s_2 ’ s_1 


+ V'{t)\\{d,d)\\B‘-^Bs-^ )-KV'{t)f q . 


From the above analysis, we deduced that 

(3.14) 2^- 2 )/ 9 «2^max(l,2- 2 «)||A 9 5|| i2 +2«( s - 1 )max(l,2-«)||A ? J|U2. 

Step 3 : the damping effect. For large enough constant K, we have 

E (Ca,(t)||(o,d)|| s ._ 2l . x5 ._ 2l ._ 1 -K2^-Vf q [t)) < 0. 
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So after simple calculations, we have 


(3.15) 


- ( ^2« (s - 2) / 9 | +c||5(i)|| J 

i?ez 


+ c^ 2 q ^ s 1 ^min(2 29 , l)max(l, 2 q )\\A q d(t)\\L 2 

<C\\F(t)\\ s .- a ,.+C\\G\\ s .- a ,.- 1 . 


Step 4 : the smoothing effect. In this step, we investigate the smoothing effect 
on d. Due to (13.151) . we just need to prove it for high frequencies only. Let k 2 = 
IIA^II 2 ,. We can easily get 

l -jk 2 q + i/||AA g J|| 2 2 =(A q G\A q d) + (AA q d\A q d) 

- (A q {u ■ Vd)\A q d) - KV'(t)k q 

Choosing q > q 0 + 1, we have ||AA g d||| 2 > & 2q \\A q d\\ 2 L2 . Hence, we get 

jk q +c2 2q k q < C'||AA g a|| i 2 + \\A q G\\ L z+CV\t)a q 2- 2 ^ s -V\\d\\B S -.,s-,. 
Summing up for q > qo + 1, we have 


- £ 29(S_1) M+ E 2^ +1 )||A 9 J|| l2 

(3.16) \9>9o + l / 9>90 + l 

<c ^ 2« s ||A 9 a|| L 2+C||G||^-2, s -i+C'W(t)||J||B S _ 2 , s _ 1 . 

<J><?0 + 1 

Taking a small positive constant e, and performing the following calculations 

(13T51) +6 (13461) . 


We will get 

A e -^( t )^(a ?rf) _ hc || S(t) || i . s+c ||J( t )| UatS+1 < C '|| J P( t )||^_ 2 , s 

+ G||G(t)|| 53 _ 2 , s _ 1 +GW(t)||J(t )||^_ 2 , s _ 1 

where 


T s (a,d) =e KV ^ ( Y^ s -Vf q +e ^ 2^ s ~^k q 

z 9>® + l 


Obviously, we know that 

J~ 5(0,6?) ~ ||o(t)||^ s _ 2iS T ||d(t)||^ s _ 2 , s _i. 

Using the above identity, further more, we can obtain 

|e-^( t )j- s (a,d)+c||a(t)||^+c||d(f)b s , 3+1 < G||F(i)|| 5 ._ a ,. 

+ C||G(f)|| ^_ 2 , a _ 1 + CV'(t)e- KV ^T s (a, d). 
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Denote C + K as C, we have 

(3 17) j t e- cv ^T s {a,d)+ce- cv ^ (||a(*)|| J& . + ||d(i)|| s .,. +1 ) 

< Ce- cv{t) (||F(t)||s s _ 2 , s + ||G(t)||B.- ai ._ 1 ) ■ 

Step 5 : the equation of fi. Cl satisfies 

d t Cl - pAQ = H - KV'(t)Cl. 

Localizing the above equation, we find 

d t A q Cl - pAA q Cl = A q H - KV\t)A q Cl. 

Hence, we can easily get 

^e~ KV ^M)\W+c2^A q mW <C\\A q H{t)\\„. 

Noting the definition of hybrid Besov space, we obtain 

As in Step 4, denote 2K as C, we have 

(3 18) +ce- cv M\\mh*,s +1 

<Ce- cv U\\H(t)\\Bs — 

Denote £ s {a,d, fl) = F s + ||H(t)||^ s _ 2 , 3 -i. Obviously we have 

£ s (a,d,£l) ~ ||a(t)||^ s _ 2 , s + ||d(t)||g s _ 2 , s -i + ||H(t)||^ 3 _ 2|3 _i. 

Summing up (13.171) and (I3.18|) . we get (13.41) . Hence, the proof is completed. □ 

Remark 3.2. For dimension N = 3, if we take s = N/2, using similar methods, we 
know that there exists a functional 

£ s (a,d,fl) ~ ||a(t)||^iv/ 2 -i,jv /2 + ||d(t)||gAr/ 2 -i + ||fl(t)||^iv/ 2 -i ) 

such that 

CL ^ e -cy(t) £ s (a,d,n)) + ce~ cvw (\\a(t)\\g N/2 + \\d(t)\\^ N/2+ i 

+ ll^Wlljv/i+i) — Ce cv ^ (\\F{t)\\^ N / 2 -i,N/2 + ||G(t), H(t)\\j 3 N/ 2 -i ), 
where C depends on p, A. 

4. Energy evolution of negative Besov norms 

In this section, we will derive the evolution of the negative Besov norms of the 
solution. The negative Besov space also used in |18j . however, we derive a different 
form of estimates for our low regularity assumption. 

Proposition 4.1. For s £ [0, N/2], we have 

\\{a{t),u(t))\\l- s < e /o(ll“llB^+l+ll«llBN/ 2 +l,^/ 2 +2 )dr|| (a0)Uo) |||_ s _ 
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Proof. From the systems we obtain 

dtAqd + divA g zt = —A g (adivu) — A q {u ■ Va), 
d t A q u — /.lAAqU — (A + /z)VdivA g u + VA g a 

= —A q (u ■ S7u ) — A q ( h(a)(fiAu + (A + ^)Vdivzz)) — A g (/(a)Va). 

Denoting (•)•) as the L 2 inner product. Applying the operator A -s to both equations 
and multiplying A~ s A q a, A -S A q u to the above two equations separately, we have 

(dtA~ s A q a\A~ s A q a) + (A _s divA g u|A~ s A g a) 

= — (A _s A g (adivzz)|A _s A g a) — (A~ S A q (u ■ Va)|A _s A g a), 

and 

((9 t A _s A g zt|A _s A g u) — fj,(A~ s AA q u\A~ s A q ii) — (A + /z)(A _s VdivA g u|A _s A g M) 
+ (A _s VA g a|A _s A g zz) = -(A~ s A q (u • Vu)\A~ s A q u) 

— (A _s A q (h(a)(fiAu + (A + /i)Vdivzi))|A~ s A g u) 

— (A~ s A q (f (a)V a)\A~ s A q u). 

Summing up the above two equalities, we find 

(41) lj t \\(A- s Aqa : A- s A q u)\\l 2 +v\\A- s VA q u\\l 2 

+ (A + /j.)||A s divA g tt ||| 2 = W\ + W 2 + W 3 + W 4 + IF 5 , 

where 

W x = — (A _s A g (adivu)|A _s A g a), W 2 = -(A~ s A q (u ■ Va)|A“ s A g a), 

W 3 = ~(A~ s A q (u ■ Vu)jA~ s A q u), W 4 = — (A _s A g (/(a)Va)|A _s A g u), 

W 5 = — (A~ s A q (h(a)(idAu + (A + /^)Vdivit))|A _s A g zt). 

Now, we give the detailed proof of the estimates about W\. 

W-l < ||A _s A g (adivM)|| i 2 ||A _,s A g a || L 2 

< C (llTadivulls-^ + HTdiv^aUg-^ + ||i?(a, divu)|| 6 -^^ Iklls-^- 
The hrst term in the bracket can be estimated by using Lemma 12.71 as follows 
||T a divw||^-s^ < C|l a lll|divu|| qn/ 2 - 

The second term in the bracket can be estimated by using Remark 12.81 as follows 
ll^divuC'll^-^ A C , ||di'vu||^o j 

< C'||divzi|| g N/ 2 1 |a|| b~ s ^- 

The third term in the bracket can be estimated by using Lemma 12.91 as follows 
||R(a,divu)|| i j-^ < C\\a\\g-^ lldivu^^. 

Combining the above three estimates with (14.21) . we find that 
(4.3) Wr < q|a|||_ JM|^ / 2 + i. 
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, mt; luiiuwmjd, 
II2 


Similar to the estimate about Wi, we can get 

W 2 < C1M|£-^||u||£-sJ|a||gAr/2 + i, W 3 < C'Hull^-^ llull^jv/2+i, 

W4 < C||a||g-^ ||u||g-s^ ||a||gAr/2+i, W5 — C\\u\\g-sJ\ a \\^-sJ\u\\^ N/2+ , 

Plugging the above estimates for W\ to W 5 to equation (ED, we will get 
^||(a(t),u(t))|||_^ + c||Vu(f)|||_^ 

: C (ll a (OII B N / 2+1 + ll M (^)ll BN/2+1,N/2+2 


< 


where s G [0, iV/2]. Finally, by Gronwall’s inequality, we obtain 

(4.4) ||(a(f),u(f))]||_ s < e /o(ll“llA~/2+i + ll“ll^/*+i^/2+2)^||( 00)Uo )||2_. . 

Hence, the proof is completed. 


□ 


Remark 4.2. From Proposition 14.11 we see that if we want to control 

■ D 2,oo 

we need to control 

(4.5) / (\\a\\g N / 2 +i \\u\\g N / 2 +i,N/ 2 + 2 ) dr. 

Jo 

When dimensional IV = 3, if the solution (p,u) belongs to the following space 

p - p G C([0, + 00 ); H 3 ), u G C([0, + 00 ); H 2 ) n i 2 (0, T; H 3 ), 

then we can control the above term (14.51) by using some basic interpolation inequal¬ 
ities. 


5. Derive optimal time decay rate 

Firstly, we give the following global well posedness results without proof for it is 
similar to [Si- 

Theorem 5.1. Suppose that N > 2. There exists a positive constant ao such that 
if po — 1 G B N l 2 ~ l ’ N / 2 + l t UQ g b n / 2 ~ 1,n / 2 yjtfh moreover 

||P 0 — l|lsN/ 2 -l,JV/ 2 +l + ||'Mo||BiV/ 2 -l,N -/2 < OO, 

then: System ED has a unique global solution ( p , u) such that 
II 1|Ic(R+,B«/2-1,N/2 + 1) + IM| C ( K+i BJV/2-1,JV/2) 

+ II P~ 1|Ili(r+,b^/2+i) + II w IIbi(r+,bn/2+1,N/2+2) < Cao- 

Applying Proposition 13. II to the system (13.21) with s = -y + 1, we will get 
(5 1) ~dt 6 Cl ^ ^ N / 2 + 1 ^ a ' d' + ce (IMIb^/w + ||(^,^)IIbn/ 2 + i ,n-/2+2) 

<Ce- cv M(\\F(t)\\z N/2 - 1 , N/2+1 + \\G(t)\\g N/2 -l,N/2 + \\H(t)\\gN/2-l,N/2) , 

where 

^Af/2+l( a , d, ~ ||a(0llB iV / 2 - 1 . iV / 2 + 1 + \W{t)\\p jN/2 -i,N/2- 
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Then, we need to give the estimates about F, G and H. We will use Lemma [2~7l 
and Lemma 12.91 frequently without mentioned. Remaindering that F = — adivit, so 
we have 

||F|| BN/2~l,N/2+l ^ ||Tadivu||gAr/2-l,AT/2 + l T II ^divu® II jjN/2-1, N/2+1 

+ \\R{a, divu)\\^ N / 2 -i,N/ 2+1 
(S ll a llijW/2-l,N/2 ||diVM||^ JV /2,N-/2 + l 
(5.2) + ||divu||^iv/2 ||a||^iv/2,Ar/2+i 

+ ||o||gN‘/2-l,N/2 ||divZt||^JV/2,N/2 + l 

< C|l a ll BN/2-1.N/2+1 ||zi||gN'/2 + l,JV/2 + 2 

< CaolMI gN/2+l,N/2+2' 


For G and FT, we give the estimates of main terms as follows 

\\u * V'u||^JV/2-l,7V/2 ||F’t i V'?i||^iV/2-l,B/2 + ||TvuW|| BN/2-l,N/2 

+ ||i?('U, Vzt) |||jJV/2-l,JV/2 

-g ~ ll U llsJV/2-l,JV/2||Vu||gJV/2 + || Vu|| gN/2 || Zt|| BN/2-1,N/2 

+ ll M ll ijlV/2—1.2V/2 || \ 7 lt|| ^2V/2 

— c|I w IIb - w / 2_i >- n / 2 II u IIb w ‘/ 2 + i i n / 2 + 2 

< c^olMI £ AT/2+1, AT/2+2 •> 


(5.4) 


||/i(a)Au|| 

fiN/2-1,N/2 <C\\a\\ BN/2 II Art || BN/2-1,N/2 

<C\\a\\ BN/2-l,N/2+l ||'w||^iV/2+l ) JV/2+2 

< CaolMI BN /2 + 1 , N/ 2+2 5 


||/(a)Va|| £jAT/2-l,AT/2 <C\\a\\ BN/2-l,N/2 | | V(l|| B N / 2 

(5.5) < C'||a||^jv/ 2 -i,jv/ 2 +i ||a||^ N / 2 +i 

< Cao||o|lBN/2+i- 

Combining estimates (15731) . (15~H) and (ESI), we get the estimates for G and H as 
follows 

(5.6) ||G||gAr/2-l,AT/2 < Ca 0 (||a||^iV/2+l + \\u\\g N /2+l,N/2+2) , 

(5.7) ||L7||^at/ 2 _i iN -/ 2 < Gao (||a||^«-/2+i + \\u\\^n/2+i,n/2+2) ■ 

Plugging (15.21) . (15.61) and EH into EU>, we find 

(5 8 ) ~dt 6 Cl (^ ^ A^ / 2 + 1 ^ a, ce CV ^ (ll a llsw/ 2 +i + \\u\\ bn/2+i,N/2+2) 

< Ca^e cv ^ (||a|| bn/2+i + ||zt|| bn/2+i,N/2+2) ■ 

Denoting £ N / 2 +\{a,d,£l) simply as £’jv/ 2 +i(^ : ), and taking a 0 small enough, we 
hnally get 

—e cv ^£n/2+i(/) + ce Cl ^ (||«|ls JV / 2 + 1 + IMIsw/2+1,N/2+2) < 0. 


(5.9) 
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of hybrid Besov space, we know that 


IMI ^W/2-l,lV/2+l ~ 

; a bn/2-1 

IMI |j1V/2-1,1V/2+1 ~ 

: M b*/ 2 - 1 

IMI BN/2+1.N/2+2 ~ 

: M B iv / 2 + 1 

By Lemma 12.61 we have 



i-i + I|m||_b jv / 2 + i > 

+ IMIgJV/2+2- 


IK*)IIb*/ 2-1 < C 'll a ( t )ll|-^ll a ( t )ll^/2+l 

||u( 0 IIbiv/ 2 -i < 


where 9\ = Af y 2 ^. 1+S ■ Similarly, we also have 

<cii W (t)iil_ L iKi)iiy/2 +1 , 

IK^IU / 2 <c||u(t)|||_^||u(t)||y; 2+2 , 

where 0 2 = jv/ 2 + 2 +s • Combining the results in Proposition 14.11 and Theorem 15.11 


we find 

H a llL“(R+,B“^) 
So we easily obtain 


+ IML 0 


(-,B“^) — ^ (jl a ollB“;: 


■ II m oIIb- s 


||a(t )||^ /2+1 > CWam^ = CWam^JJli^ 1 , 

( 5 - 10 ) \\u(t)\\z N/2+1 > CMtn^-r = C\\u{t)\^^~\ 

> c\\um^f 2 = c\\u{t)\^j^. 

By taking ao small enough in Theorem 15.11 we can assume 

ll a W II BN/2-1,N/2+1 ^ lj ||n(t) HjjiV/2-l,iV/2 ^ 1- 

2 , we find 


Notice that 1 + ^4^ > 1 + 

Plugging (15.111) into (15.91) . we obtain 
—e cv (t ^Ar/ 2 +i(f) 

+ c(e ^ (||a(f)||^jv/ 2 -i,N/ 2 +i 


\Ht)\\g N/2+: 

ll U (^)ll BN/2+2 


i>GK*)||^f-, 


^ ' H M (£)llijiV/ 2 -l,N/ 2 )^ S < 0 . 

Since £N/ 2 +i(t) ~ l|a(f)||^^/ 2 -i,jv/ 2 +i + \\u(t)\\^ N/2 _ 1 , N/2 , we finally get 

j t e~ cv ^e N/2+1 (t) + c(e- cl ' (, >^ /2+ i(t)) I+TOfcT < 0 . 


(5.12) 


—o ^A/2+iv^y i ° wv/ 2 +: 

Solving this differential inequality, we could obtain 

£n/2+i (t) < Ce cv ^ ( &N/2+ 


( 0 ) «/ 2 + s -! + 


2 C 


iV/2 + s-l 


iV/2 + . 


- 1 
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From Theorem ED we know that V(t) is bounded by the initial data, there exists 
a constant C such that 

, N/2-l+s 

(5.13) ||Q-(t) ||^AT/2-l,iV/2+l + ||lA(t)||_gJV/2-l,JV/2 ^ C(\ + t) 2 

Due to the relationship between Sobolev space and Besov space, we obtain 

||A iV / 2 - 1 a(t)|| i2 + ||A w / 2 - 1 u(t)|| i2 
<C (||a(t)||g iV / 2 -i + ||u(t)||^ iV /2-i) 

<C (\\a(t)\\g N / 2 -l,N/ 2 +l + \\u(t)\\g N / 2 -l,N/ 2 ) 

N/2-l+s 

<(7(1 + 1) ——. 

For £ £ (— s, N/2 — 1), we have 

IK*)IIb< <ll«Wll|- L ll a Wliy / 2 -i, 

where 9 = (15.141) . we then obtain 

N/2- 1 — l N/ 2-1-1 

IK*)IIb< < llaWH^IIaWII ^^ 5 
< C(l + t)~*¥. 

Similar to the above analysis for aft), we also have the following estimates for u(t) 

ll u (£)ll b 1 — ^(1 + t) 2 ■ 

At this stage, we complete the proof of Theorem ll.il 

Remark 5.2. In order to get the Theorem 11.61 we need the following theorem which 
is proved in m- 

Theorem 5.3. Suppose that N > 2. There exists a positive constant ao such that if 

Pq — 1 £ B N / 2 ~ 1,N / 2 , Uq £ B^/ 2 -l wmoreover 

\\pO ~ 1||biv/ 2 -i,iv/ 2 + ||uo|IbJV/ 2 -i ^ «0) 
then: System has a unique global solution (p, u ) such that 

\\P ~ l|lc(R+,B JV / 2 - : L JV /2) "I” IMIcyR+.B-Wa-i) 

+ IIp~ 1|Ili(r+,b n / 2 ) + IMIl 1 (r+,b a, / 2 + 1 ) — Ca 0 - 

With this theorem, Remark l3.21 Remark 14.21 and the result in paper [2Tj, we can 
mimic the procedure of the last section to get our Theorem 1 1.61 
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